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ABSTRACT 

Observations of ionized and neutral gas outflows in radio-galaxies (RG) suggest that AGN radio jet feedback 
has a galaxy-scale impact on the host ISM, but it is still unclear how the molecular gas is affected. Thus it 
is crucial to determine the physical conditions of the molecular gas in powerful RG to understand how radio 
sources may regulate the star formation in their host galaxies. We present deep Spitzer IRS high-resolution 
spectroscopy of 8 nearby RG that show fast Hi outflows. Strikingly, all of these Hi-outflow RG have bright 
H2 mid-IR lines that cannot be accounted for by UV or X-ray heating. This strongly suggests that the radio 
jet, which drives the Hi outflow, is also responsible for the shock-excitation of the warm H2 gas. In addition, 
the warm H2 gas does not share the kinematics of the ionized/neutral gas. The mid-IR ionized gas lines (with 
FWHM up to 1250 km s" 1 for [Nell] 12. 8 /im) are systematically broader than the H2 lines, which are resolved 
by the IRS in tn 60% of the detected lines (with FWHM up to 900 km s _1 ). In 5 sources, 3C 236, 3C 293, 
3C 459, 4C 12.50 and PKS 1549-79, the [Nell] A 12. 8 /xm line, and to a lesser extent the [Nelll]A 15.5/an 
and [Nev]A 14.3/xm lines, clearly exhibit blue-shifted wings (up to -900 km s" 1 with respect to the systemic 
velocity) that match well the kinematics of the outflowing Hi or ionized gas. The H2 lines do not show these 
broad wings, except tentative detections in 4C 12.50, 3C 459 and PKS 1549-79. This shows that, contrary to the 
Hi gas, the H2 gas is inefficiently coupled to the AGN jet-driven outflow of ionized gas. While the dissipation 
of a small fraction (< 10%) of the jet kinetic power can explain the turbulent heating of the molecular gas, our 
data show that the bulk of the warm molecular gas is not expelled from these galaxies. 

Subject headings: galaxies: jets, galaxies: kinematics and dynamics, galaxies: evolution, galaxies: ISM, shock 
waves, ISM: general, ISM: jets and outflows, molecular processes, turbulence 
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1. INTRODUCTION 

Active galactic nucleus (AGN) feedback is recognised to 
have an important effect on galaxy evolution. Widely intro- 
duced in numerical simulations of galaxy evolution to clear 
the circum-nuclear gas and halt the growth of super-massive 
black holes (Silk & Rees 1998; Di Matteo et al. 2005), this 
mechanism would explain the correlations found between the 
BH mass and e.g. the bulge mass (e.g. Ferrarese & Merritt 
2000), and prevent the formation of too many massive galax- 
ies in the early universe (e.g. Thomas et al. 2005). 

The co-evolution of massive black holes and their 
host galaxies has been well established observationally 
in samples at various redshifts (e.g. Tremaine et al. 2002; 
Alexander et al. 2005). Nevertheless, detailed observations of 
individual active galaxies in the nearby universe, where phys- 
ical processes in the central region can be studied in detail, are 
essential for answering key questions about the role of AGN 
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in galaxy evolution, such as what is the magnitude of AGN 
feedback (e.g. gas outflow rates) and what is the main driving 
mechanism of AGN feedback (e.g. quasar wind, radio jets, 
circum-nuclear starbursts). 

Most of the existing evidence for AGN radio jet feed- 
back on the scales of galaxy bulges comes from obser- 
vations of outflows of neutral gas (Morganti et al. 2003, 
2005b; Emonts et al. 2005; Lehnertet al. 2011) and ionized 
gas (Nesvadba et al. 2006, 2008; Holt et al. 2008) in radio 
galaxies. These observations suggest that radio jets are an 
efficient mechanism to convert the energy output of the AGN 
into an energy input into the ISM, and that radio sources may 
regulate gas cooling in early -type galaxies (Best et al. 2005, 
2006; Donoso et al. 2009). 

One of major open questions about estimating mass out- 
flow rates is whether the ionized and Hi gas traces the domi- 
nant phase of the wind. Therefore it is crucial to compare the 
masses and kinematics of each gas phase. It has only very re- 
cently been recognized that starburst and AGN-driven winds 
also include molecular components (e.g. Veilleux et al. 2009; 
Feruglio et al. 2010; Alatalo et al. 201 1). Fischer et al. (2010) 
and Sturm et al. (2011) reported the discovery of molecular 
outflows in nearby ULIRGs (e.g. Mrk 231) through the de- 
tection of a series of OH and H2O lines seen in absorption 
against the bright IR dust continuum with Herschel/PACS in 
the A78 - 79 /im and A 1 19 - 121 /im ranges. However, esti- 
mates of outflow rates from absorption lines suffer from major 
uncertainties regarding e.g. the geometry, and the degeneracy 
between covering fractions and column densities. Moreover, 
most of these sources are composite: the co-existence of a 
starburst, quasar, and radio source makes it difficult to infer 
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what is driving the wind. Starbursts (Rupke et al. 2005), radi- 
ation pressure from the quasar (Feruglio et al. 2010), and even 
mechanical interactions with the radio source (Reynolds et al. 
2009) have all been proposed. 

Spitzer IRS spectroscopy opened a new perspective from 
which to study the impact of the injection of kinetic energy 
on the formation of molecular gas and regulation of star for- 
mation in nearby radio galaxies. Ogle et al. (2010) found that 
30% of a sample of 55 nearby 3C radio galaxies have unusu- 
ally bright mid-IR line emission from warm (10 2-3 K) H2 gas, 
with weak tracers of star formation (SF) (e.g. PAHs). The 
H2 to PAH luminosity ratio is more than 10 times larger than 
what is expected from UV and X-ray photon heating. We pro- 
pose that the H2 luminosity is associated with the dissipation 
of a small fraction of the mechanical energy of the radio jet 
(Nesvadba et al. 2010; Ogle et al. 2010). These H 2 -luminous 
galaxies lie off the Schmidt- Kennicutt relationship, indicating 
that star-formation may be suppressed relative to nearby nor- 
mal star-forming galaxies (Nesvadba et al. 2010). Numerical 
simulations also suggest that radio jets appear to be efficient at 
injecting kinetic energy into the ISM (e.g. Wagner & Bicknell 
201 1), possibly inhibiting star formation. 

This paper reports on the mid-IR Spitzer IRS (Houck et al. 
2004) spectroscopy of 8 nearby powerful radio galaxies in 
which fast (« 1000 km s" 1 ) outflows of ionized and neutral 
gas have been detected (Morganti et al. 2005b). These galax- 
ies are ideal targets to study the impact of AGN feedback 
on the different phases of the ISM, and in particular on the 
molecular gas, and thus on star-formation. They present an 
environment where winds can clearly be associated with the 
mechanism driving them, namely the radio jet. How is the 
kinetic energy of the outflow dissipated? Is the molecular 
gas entrained in the ionized/neutral outflow? If yes, what are 
the physical processes that control the dynamical coupling be- 
tween the tenuous gas outflow and the dense molecular gas? 

Sect. 2 and 3 present the galaxy sample and details on the 
data reduction and spectral analysis. Then we discuss spec- 
troscopic diagnostics (sect. 4), focusing on the H2 gas, and 
we compare the kinematics of the ionized, Hi, and H2 gas in 
sect. 5. Sect. 6 discusses the excitation mechanisms and origin 
of the H2 gas, and proposes an interpretation of the observed 
gas heating and kinematics of the multiphase ISM gas of the 
host galaxy. 

2. THE SAMPLE OF HI-OUTFLOW RADIO GALAXIES 

Except PKS 1549-79, all of the observed sources are part 
of the Morganti et al. (2005b) sample. Most of these sources 
are bright compact radio sources, with sizes ranging from 
0.01 to 40 kpc. These Hi-outflow objects appear to be- 
long to the minority of powerful radio galaxies (w 35%) 
that show evidence for recent star formation at optical wave- 
lengths (Tadhunter et al. 2011, based on spectral synthesis 
modelling) and/or mid-IR wavelengths (Dicken et al. 2011, 
based on PAH detection). Most such objects are compact ra- 
dio sources (e.g. CSS/GPS/CSO) or have unusually strong 
compact steep spectrum cores (Tadhunter et al. 2011). The 
star formation in these objects is probably linked to the pres- 
ence of an unusually rich ISM in the nuclear regions of the 
host galaxies, perhaps accreted in gas-rich mergers. In such 
objects the radio jets are expected to interact particularly 
strongly with the dense circumnuclear gas. Therefore these 
are just the type of objects in which we might expect the jets 
to have a direct impact on the molecular gas. 

Some properties of the sample are listed in Table 1. The 



stellar masses are derived from the K-band luminosities ac- 
cording to Marconi & Hunt (2003): 

log 10 (M*) =-2.3 + 1.21 xlog m (L K ) , (1) 

where M± is the bulge stellar mass in M Q and Lk the TT-band 
luminosity in L Q . The unabsorbed 2-10 keV AGN X-ray 
fluxes are taken from the literature, except 3C 236, 3C 459 and 
PKS 1549-79 which were unpublished. For 3C 236, we re- 
processed and extracted the Chandra ACIS spectrum. We find 
that the spectrum is well fitted with a single power law of pho- 
ton index 1 .67, absorbed by a column of A^h = 2.3 x 10 23 cm" 2 , 
consistent with the presence of a dust lane close to the nu- 
cleus (deKoffetal. 2000). The 3C 459 and PKS 1549-79 
X-ray fluxes were provided by P. O'Brien and M. Hardcastle 
(private communications). The jet kinetic power is estimated 
from the 178 MHz flux density, according to the formula 10 
derived in Punsly (2005). 

These galaxies all show very broad (up to 1000 km s" 1 ) Hi 
absorption profiles, with highly blue-shifted Hi gas with re- 
spect to their optical systemic velocities, which indicates a 
fast outflow of neutral gas (e.g. Morganti et al. 2003). The 
associated Hi mass outflow rates are up to ps 60M q yr , 
up to 2 orders of magnitude higher than the mass outflow 
of ionized gas (see e.g. Emonts et al. 2005). These outflow 
rates are estimated from the gas column density, the out- 
flow velocity, and the radius within which the flow originates 
(see e.g. Heckman 2002, for a review). In some cases, 
like 3C 305 (Morganti et al. 2005a), IC 5063 (Oosterloo et al. 
2000; Morganti et al. 2007) and 3C 293 (Morganti et al. 2003; 
Emonts et al. 2005), this radius is estimated from high- 
resolution radio and optical observations, where the outflow is 
spatially resolved against the radio continuum, and co-spatial 
with regions where the radio emission in the propagating ra- 
dio jet is enhanced. For the other sources, the radii are more 
uncertain. It is assumed that the outflow is coming from the 
brightest radio continuum region, which in principle provides 
only a lower limit since the Hi absorption can only be traced 
in regions where the continuum is bright. We checked that 
these radii are consistent with optical spectroscopy, but note 
that these radii estimates only provide an order of magnitude. 
Some of the main radio and Hi properties are summarized in 
Table 2. In sect. 5.4, we discuss in more detail the Hi kine- 
matics, and compare them to the kinematics of the molecular 
and ionized gas phases. 

3. SPITZER OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

The IRS has two high-resolution modules, short high" 
(SH) and long-high 12 (LH), and four low-resolution modules, 
short-low (SL1 & SL2) and long-low (LL1 & LL2), spanning 
5.2-38 fim. In low resolution mode (A/AA » 57- 127), we 
used publicly available archival data for 7 of the 8 sources. 
Only PKS 1549-79 was not previously observed at low res- 
olution, hence we observed this source as part of our Spitzer 
Cycle 4 program p40453. At high resolution (A/AA « 600), 
we observed 5 sources (3C 236, 3C 305, 3C 459, PKS 1549- 
79, and IC 5063) as part of our program p40453, and we used 
archival data for the remaining 3 sources (3C 293, 4C 12.50, 
and OQ 208) (p30877). The exposure times for each mode are 

10 This formula uses the 151 MHz flux density. When a direct measure- 
ment at 150 MHz was not available, we derived it from the 178 MHz measure- 
ment (taken from NED) by fitting the measured radio SED with a power-law. 

11 4.7" X 11.3" slit, observed wavelengths from 9.9 to 19.6/^m 

12 11.1" X 22.3" slit, observed wavelengths from 18.7 to 37.2 ^m 
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Table 1 

Properties of Hi-outflow radio-galaxies observed with Spitzer IRS 



Object 


RA a 


Dec 1 * 


b 


D L C 


K-mag 


M* d 


M(H 2 f 


M(H 2 ) 


i-24/im f 


logL x g 


logZ. x 


L(H 2 )/L x h 




J2000.0 


J2000.0 




[Mpc] 




[M Q ] 


[M ] 


ref. 


[Lq] 


[erg s~'] 


ref. 




3C 236 


10 06 01.7 


+34 54 10.4 


0.1004 


449 


12.25 


7.4E10 


< 8E9 


(1) 


1.94±0.15 E10 


43.0 




0.054 


3C 293 


13 52 17.8 


+31 26 46.5 


0.0448 


195 


10.84 


4.6E10 


1.5E10 


(2) 


6.26 ±0.75 E09 


42.8 


(4) 


0.063 


3C 305 


14 49 21.8 


+63 16 15.3 


0.0418 


178 


10.64 


4.8E10 


2.0E9 


(1) 


7. 13 ±0.45 E09 


41.2 


(5) 


1.78 


3C 459 


23 16 35.2 


+04 05 18.1 


0.2199 


1046 










5.36±0.32E11 


43.2 




0.13 


4C +12.50 


13 47 33.4 


+12 17 24.2 


0.1234 


551 


12.21 


1.3E11 


1.5E10 


(3) 


1.20±0.04E12 


43.3 


(6) 


0.16 


IC 5063 lobe 


20 52 02.1 


-57 04 06.6 


0.0113 


45.3 


8.75 


1.6E10 






3.74 ± 0.08 E10 


42.9 


(7) 


0.008 


OQ 208 


14 07 00.4 


+28 27 14.7 


0.0766 


337 


11.52 


5.2E09 


1.4E10 


(1) 


2.31 ±0.03 Ell 


42.7 


(8) 


0.08 


PKS 1549-79 


15 56 58.9 


-79 14 04.3 


0.1521 


699 










1.31 ±0.04 E12 


44.7 




0.005 



References. — (1) Ocana Flaquer et al. (2010); (2) Evans et al. (1999); (3) Evans et al. (2002); (4) Ogle et al. (2010); (5) Massaro et al. (2009); (6) O'Dea et al. 
(2000); (7) Lutz et al. (2004); (8) Guainazzi et al. (2004) 

a Positions of the sources used for the Spitzer IRS observations. Except for 3C 305, IC 5063 and 4C 12.50, the positions of the radio nuclei are taken from NED 
(see Sect. 3 for details). 

b Redshift derived from optical line measurements (see references for the systemic velocities in Table 2). 
c Luminosity distance assuming Ho = 73km s Mpc~' and f2,„ = 0.27. 
d Stellar mass derived from the K-band luminosity. 

c H2 gas masses estimated from CO(1-0) measurements. For 4C 12.50, only the H2 mass in the western lobe is quoted. 

f Rest 24 fjxa luminosity measured with Spitzer IRS. We averaged the flux over 23 - 25 fim, after removal of the [NeV] 24.3 1 pm when present. 
g Unabsorbed 2—10 keV nuclear X-ray luminosity in erg s ■ For 3C 236 for which we extracted the X-ray spectrum from the Chandra archive (see sect. 6.1). 
The PKS 1549-79 X-ray flux is from Chandra data (Paul O'Brien, private communication) and agrees with the Suzaku value (4.5 X 10 44 erg s" 1 ). The 3C 459 
X-ray flux is from XMM data (Martin Hardcastle, private communication). 

h Ratio of the H2 line luminosity (summed over the S(0) to S(3) rotational transitions) to the unabsorbed 2—10 keV nuclear X-ray luminosity. 



Table 2 

Radio and Hi properties of outflow radio-galaxies observed with Spitzer IRS 



Object 


size a 


size 


^178MHz 


Qje, h 


T 


A'(Hl) 


M(Hl) c 


v d 


I'sys 


Vout 


M 


radius 




[kpc] 


ref. 


[Jy] 


[ergs '] 




[10 21 cnr 2 ] 


[M ] 


[kms- 1 ] 


ref. 


[kms-'j 


[Meyr 1 ] 


[kpc] 


3C 236 


6.0 


(1) 


15.7 


7.6 x 10 45 


0.0033 


5.0 


1.2E8 


30129 


(9) 


750 


47.0 


0.5 


3C 293 


5.0 


(2,4) 


17.1 


5.1 x 10 45 


0.0038 


6.0 


1.0E8 


13450 


(2) 


500 


56.0 


1.0 


3C 305 


10.0 


(3) 


13.8 


6.2 x 10 45 


0.0023 


2.0 


4.7E7 


12550 


(7) 


250 


12.0 


1.0 


3C 459 


39.0 


(5) 


30.8 


2.5 x 10 46 


0.0005 


0.75 


1.8E7 


65990 


(10) 


300 


5.5 


1.0 


4C 12.50 


0.15 


(6) 


4.6 


2.5 x 10 45 


0.0017 


2.6 


7.0E6 


37027 


(10) 


600 


21.0 


0.02-0.2 


IC 5063 


8.0 


(7) 


6.3 


2.0 x 10 45 


0.0120 


10.0 


2.4E8 


3400 


(11) 


350 


35.0 


0.4 


OQ208 


0.01 


(6) 


0.12 


5.2 x 10 43 


0.0057 


8.3 


2.0E4 


22985 


(12) 


600 


1.2 


0.01 


PKS 1549-79 


0.4 


(8) 


12.6 


7.7 x 10 45 


0.02 


40.0 


1.0E8 


45628 


(10) 


250 


30.0 


0.1-1.0 



References. — (1) deKoffetal. (2000); (2) Emonts et al. (2005); (3) Morganti et al. (2005a); (4) Zirbel & Baum (1998); (5) Nilsson et al. (1993); 
(6) Xiangetal. (2002); (6) Heisler & Vader (1994); (7) Morganti et al. (2005b); (8) Holt et al. (2006); (9) Hill et al. (1996); (10) Holt et al. (2008); (11) 
Morganti et al. (1998); (12) Marziani et al. (1993) 

a Size of the radio source. Note that for 3C 293, this is the size of the inner radio source. 
b Jet cavity kinetic luminosity estimated from the 178 MHz flux (see text, Sect. 2). 
c Mass of Hi gas in the outflowing component (within the radius quoted in the table). 
d Systemic velocity of the galaxy (see text for details). 

c Outflow velocity, estimated as half of the full width at zero intensity of the blue-shifted part of the broad Hi component. 
f Radius of the Hi outflow region. Except for PKS 1549-79, the radii values are from Morganti et al. (2005b). 

listed in Table 3, and were chosen to reach a signal-to-noise 
ratio of at least 10 for the H2 S(l) 17 /im line. 

Even if the outflow of neutral gas is occurring at a kpc-scale 
distance from the nucleus, this is still spatially unresolved 
with the IRS in high-resolution mode. However, two sources 
(3C 305 and IC 5063) exhibit a more extended inner radio 
structure, and we placed the slit slightly off-nucleus, coinci- 
dent with a bright peak in the radio continuum, at a distance 
of 2 arcsec from the core. Note that the core is still included 
in all of the IRS slits. The slit positions are given in Figure 1. 
4C 12.50 is a system of two interacting galaxies (Axon et al. 
2000), and the SH slit was positioned only to include the west- 
ern companion. For the other sources, the slits were centred 
on the nucleus. 



Table 3 

Spitzer IRS observation log 



Object 


SL 




LL 




SH 


LH 


3C 236 


13 X 14 


5 


X 


14 


4 


X 


120 


10x60 


3C 305 


13 x 14 


5 


X 


14 


4 


X 


120 


10x60 


3C 293 


3 x 120 


3 


X 


120 


3 


X 


120 


10x60 


3C 459 


16 x 14 


2 


x 


30 


4 


X 


120 


10x60 


4C +12.50 


3 X 14 


2 


X 


30 


6 


X 


30 


4x60 


OQ 208 


2x60 


2 


X 


14 


4 


X 


120 


3x60 


IC 5063 lobe 


2x 14 


1 


X 


30 




X 


30 


4x60 


PKS 1549-79 


2x60 


3 


X 


30 


4 


X 


120 


5 x 60 



Note. — Exposure times for the IRS modules given as 
number of cycles X ramp duration in seconds. 



3.1. Spitzer data reduction 



The galaxies were observed in staring mode, and placed at 
two nod positions along the slit (1/3 and 2/3 of the length 
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IC5063 

HST WFC2 [0111] 5007A image 




Figure 1. Footprints of the high-resolution IRS slits (SH and LH) for the 3C305 and IC5063 galaxies. For these two sources, the radio observations resolve the 
Hi outflow against the radio jet, at a distance of about 2 arcsec from the nucleus. Therefore, the slits were centred slightly off-nucleus, at the location of the 
broad, blue-shifted Hi absorption associated with the outflow (indicated by the orange circle). The position of the radio core is marked with a green cross. 



in IRS staring mode, and the SMART package 16 was used to 
perform optimal extraction. We checked on 3C 293 that the 
SPICE and SMART extractions give very similar results. 

3.2. Spectral analysis 

To measure the fluxes and linewidths of emission lines, we 
first fitted the spectra with the IDL PAHFIT(Smith et al. 2007) 
tool, which decomposes the spectra into starlight, thermal 
blackbody emission (extinguished by a uniform dust screen 
if specified), resolved PAH features (fitted by Drude profiles), 
and emission lines. Originally designed to fit low-resolution 
Spitzer IRS spectra, we adapted the routine to handle high- 
resolution data. In particular, we allow the emission lines to 
be spectrally resolved and we fit them by Gaussian profiles. 

We use the PAHFIT decompositions to measure the fluxes 
of the PAH complexes, and to produce starlight-free and dust- 
free spectra. Except for OQ 208 and PKS 1549-79, we do not 
include any extinction in the fit. For OQ 208, 4C 12.50 and 
PKS 1549-79, we introduced the capability of fitting silicate 
absorption features at 9.7 and 18 /iin. Then we ran PAHFIT a 
second time to fit accurately the emission lines. We also fit- 
ted the lines with Gaussian profiles using the non-linear least- 
squares IDL fitting routine MPFIT (Markwardt 2009) on the 
continuum-free spectra. The two methods agree very well and 
give similar line fluxes (within the uncertainties). We find 
significant differences (up to a factor of 3 in line flux) when 
trying to measure the line fluxes in one go with PAHFIT, the 
many free parameters making the fit inaccurate. 

For each of undetected emission lines, we computed the 
flux of an unresolved Gaussian line with a peak flux equal to 
twice the rms noise in a 1 ^in-wide band centred at the wave- 
length of the line. We took this flux value as a 2 a upper limit. 

We also compute the 24 pm rest-frame fluxes by averaging 
the IRS line-free spectrum over a narrow band (23-25 pm). 

13 http://irsa.ipac.caltech.edu/data/SPlTZER/docs/dataanalysistools/tools/cateStt|^rrowband 24 pm fluxes are listed in Table 1. 

14 http://irsa.ipac.caltech.edu/data/SPlTZER/docs/dataanalysistools/tools/irsclean/ 

15 http://irsa.ipac.caltech.edu/data/SPlTZER/docs/dataanalysistools/tools/cuDfehlt'p://isc. astro.cornell.edu/IRS/SmartRelease 



of the slit). The data were processed by the Spitzer Science 
Center S18.7 pipeline. We started the data reduction from 
the basic calibrated data products, which include corrections 
for flat-fielding, stray light and non-linearity of the pixels re- 
sponse. All but one (4C 12.50) galaxies had dedicated sky 
background (off) observations for all high-resolution mod- 
ules. When available, the 2D off spectra were subtracted from 
the on observation for each nod position. For 4C 12.50, the 
LH module was used to remove the background emission. For 
the LH observations, we corrected for the time- and position- 
dependence of the dark current using the DARK_SETTLE 
routine 13 . 

The cleaning of rogue pixels was performed using the 
IRSCLEAN routine 14 . Rogue pixel masks of the correspond- 
ing IRS campaigns were used, and a further manual clean- 
ing was done for each nod position and module. The Spitzer 
IRS Custom Extractor (SPICE) software was used to extract 
the one-dimensional SH and LH spectra, with optimal extrac- 
tion. The individual orders within the high-resolution mod- 
ules were very well aligned in flux and no scaling was applied 
between the orders. Then the orders were trimmed (by re- 
moving the end of each order, where the response drops) and 
stitched together. Finally a median combination was used to 
produce the final stacked spectra. 

For all but one source (4C 12.50), the continuum levels of 
the SH and LH matched and no scaling was necessary. For 
4C 12.50, we scaled the SH spectrum by a factor of 1.1 to 
match the LH module, since it is the most accurate due to its 
background subtraction. 

The 3C236 and 3C305 low-resolution spectra were ex- 
tracted from spectral maps with the CUBISM 15 software. For 
the other galaxies, the low-resolution observations were done 
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4. SPECTROSCOPIC DIAGNOSTICS 

The full high-resolution (SH+LH, i.e. observed wave- 
lengths 9.9- 37.2 yum) spectra are shown in Figure 2, and the 
low -resolution spectra are presented on Figure 3. The loca- 
tion of the detected PAH complexes and emission lines are 
marked on each spectrum. Overall, the spectra are of very 
good quality, with noise levels expected for such exposure 
times. Some minor spurious features still remain, in partic- 
ular in the LH spectra (noisier than the SH), which are due 
to bad pixels or instrumental artefacts that were not cleaned 
during the data reduction. In appendix, we provide zooms on 
the SH wavelengths (Figure 9), and on all individual detected 
emission lines (Figure 10 and 11). 

4.1. Molecular hydrogen 

Powerful emission lines from H2 are detected in all of the 8 
high-resolution spectra. The H 2 0-0 S(l) 17/xm and 0-0 S(3) 
9.7 /j,m lines are detected in all radio-galaxies. The pure ro- 
tational H2 line fluxes (or their 2 a upper limits) are listed in 
Table 4. In some cases, like 3C 236, 3C 293, and4C 12.50, the 
H2 S(l) line dominates the SH spectrum (Figure 9). These H2 
strengths are comparable to those observed in the Ogle et al. 
(2010) sample of Ha-bright radio-galaxies, in some galaxy 
collisions, like Stephan's Quintet (Cluver et al. 2010) and the 
Taffy galaxies (Peterson et al. submitted). 

To compute the masses of the warm (> 100 K) H2, we fol- 
low two different approaches. The first one is similar to that 
described by Ogle et al. (2010) and consists of assuming a 
thermal distribution of the Ht levels. We constructed excita- 
tion diagrams, shown on Figure 12, by plotting the logarithm 
of the column densities of the upper H2 levels divided by their 
statistical weights, ln(N v j/gj), against their excitation ener- 
gies, E v j/ks, expressed in K. For a uniform excitation tem- 
perature, the values ln(N v j/gj) should fall on a straight line 
plotted versus E v j/ks, with a slope proportional to !T e ~'. In 
a situation of local thermal equilibrium (LTE), the excitation 
temperature T exc equals that of the gas. We fitted the data with 
two or three excitation temperature components, although the 
gas is at a continuous range of temperatures. We constrained 
the temperature range to be 100 < T < 1500 K. At each fitting 
iteration, the statistical weights of the ortho transitions were 
adjusted to match the LTE value. 

The results are given in Table 5, where we list the fitted H2 
excitation temperatures, ortho-to-para ratios, column densi- 
ties, masses, and total luminosities for each temperature com- 
ponents. The bulk of the warm H2 is constrained by the lowest 
temperature component, i.e. by the S(0) and S(l) lines. Thus, 
if the S(0) line is not detected, we include its upper limit in the 
fit as a 2 a detection (with 1 a error bar) and quote the model 
results as upper limits. In all galaxies where the S(0) line is 
detected, we measure warm H2 masses, from « 3 x 10 8 M Q 
in 3C 236, up to « 3 x 10 10 M Q in 4C 12.50. Similar ex- 
tremely large warm H2 masses were detected in 3C 433 and 
3C 436 (Ogle et al. 2010). The H2 column densities range 
from < 4 x 10 20 to « 2 x 10 22 cm -. These column densities 
depend on the assumed size of the emission region. 

The second approach to fit the H2 line fluxes and derive 
the H2 physical parameters is similar to that described in 
Guillard et al. (2009) and Nesvadba et al. (2010). It assumes 
that the H2 emission is powered by the dissipation of mechan- 
ical energy in the molecular gas (see sect. 6. 1 for a discussion 
of the H2 excitation mechanisms). We model this dissipa- 
tion with magnetic shocks, using the MHD code described 



in Flower & Pineau Des Forets (2010). The gas is heated to a 
range of post-shock temperatures that depend on the shock ve- 
locity, the pre-shock density, and the intensity of the magnetic 
field (which is perpendicular to the shock propagation). We 
use a grid of shock models (varying shock speeds) similar to 
that described in Guillard et al. (2009), at pre-shock densities 
«h = 10 3 cm" 3 and 10 4 cm" 3 . The initial ortho-to-para ratio 
is set to 3, and the intensity of the pre-shock magnetic field is 
30 y«G. The H2 line fluxes are computed when the post-shock 
gas has cooled down to a temperature of 100 K. At a given 
pre-shock density, the shock velocity is the only parameter 
we allow to vary. 

A combination of two shock velocities is required to match 
the observed H2 line fluxes. We show the best-fitting shock 
combination in Figure 13 and 14 for the two pre-shock densi- 
ties. These fits are not unique, but they provide an estimate of 
the range of shock velocities and pre-shock densities needed 
to reproduce the H2 excitation. This phase-space is well con- 
strained when six or more H2 lines are detected. The lowest 
density and shock velocity is required to fit the low-excitation 
lines S(0) and S(l), which probe the bulk of the warm Hz 
mass, whereas the high density and shock velocity is needed 
to fit the high-excitation lines, above S(4). In reality, the H2 
emission arises from a distribution of densities and shock ve- 
locities. With the present data we cannot exclude the presence 
of very high density gas («h > 10 5 cm" 3 ). 

As discussed in Guillard et al. (2009), the H2 masses are 
derived by multiplying the gas cooling time (down to 100 K) 
by the gas mass flow (the mass of gas swept by the shock per 
unit time) required to match the H2 line fluxes. The shock 
model parameters, gas cooling times, mass flows, and warm 
H2 masses are quoted in Tables 6 and 7. The warm H2 masses 
derived from the shock modeling are larger than from the LTE 
model. This is mainly because of the different values of the 
H2 ortho-to-para ratio, and because at «h = 10 3 cm" 3 , the S(0) 
and S(l) lines are not fully thermalized. 

4.2. Thermal dust continuum and aromatic features (PAHs) 

Our sample reveal two types of spectra, those with a flat 
thermal dust continuum (3C 236, 3C 293, 3C 305), and those 
with a steep rising continuum at long wavelengths (3C 459, 
4C 12.50, IC 5063, OQ 208, and PKS 1549-79). There 
is a remarkable similarity between the spectra of 3C 236, 
3C 293 and 3C 305 with the 3C 326 radio-galaxy (Ogle et al. 
2007), the Stephan's Quintet (Cluver et al. 2010) and the 
Taffy galaxies (Peterson et al., submitted), where the dust 
continuum is very weak compared to emission lines. The 
narrow -band 24 /im luminosities (see sect. 3.2 and Table 1), 
measured on the LH spectra, span three orders of magnitude, 
from 5.4 x 10 9 L (3C 305) to ULIRG-like luminosities of 
1.1 x 10 12 L (4C 12.50). For 4C 12.50 and OQ 208, our 
measurements are in agreement with those of Willett et al. 
(2010) on low-resolution IRS data. 

Emission from aromatic infrared bands attributed to Poly- 
cyclic Aromatic Hydrocarbons (PAHs) excited by UV pho- 
tons are detected in all galaxies. In general, in the high- 
resolution data, the 11.3 /.tin band is the most prominent fea- 
ture. The fluxes of the PAH complexes, given in Table 8, are 
measured with PAHFIT on the emission line-free spectra (see 
sect. 3.2), by combining the fluxes of the blended features that 
contribute to the main complexes (around 8, 11.3, 12.6 and 
17 fim). 

All of the radio-galaxies have a 11.3//m PAH to 24 /im 



6 



P. Guillard et al. 




4 ^3C293 (z = 0.04483) 



15 20 

Rest wavelength [fin] 



0.1 r„ 




15 20 25 

Rest wavelength [fin] 



, _ 3C305 (z = 0.04183) 



3 0.6 - 



0.2 -ST - L W2 



^t.iiiiiuiij 




15 20 25 

Rest wavelength [fin] 



3 . 3C459 (z = 0. 21997) 




15 20 

Rest wavelength [fin] 














_ IC5063 


lobe (z=0.01 133) 


c 












>)e V] 








> 












z 



































20 25 30 

Rest wavelength [fin] 



20 25 
Rest wavelength [fin] 




15 20 25 

Rest wavelength [fin] 



Figure 2. Spitzer IRS high-resolution (SH+LH) spectra of the Hi -outflow radio-galaxies. The locations of PAH features and emission lines are marked for each 
object. 
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Figure 3. Spitzer IRS low-resolution (SL+LL) spectra of the Hi-outflow radio-galaxies. The locations of PAH features and emission lines are marked for each 
object. The low-resolution spectra for 3C 293, OQ 208 and PKS 1549-79 were already presented by Ogle et al. (2010) and Willett et al. (2010). 
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Table 4 

H2 line fluxes from Hi-outflow radio galaxies 



Source 


H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 


H 2 S(4) 


H 2 S(5) 




S(6) 


H 2 S(7) 




28.22 fan 


17.03 Aim 


12.28 A*m 


9.66 fim 


8.03 /xm 


6.91 fj,m 


6.1 


1 


5.51 


3C 236 


0.20 (0.03) 


0.95 (0.07) 


0.32 (0.05) 


0.64 (0.05) 


< 1.04 


< 1.10 


< 


0.23 


<0.31 


3C 293 


0.79 (0.06) 


3.46 (0.29) 


1.61 (0.10) 


2.50 (0.16) 


1.35 (0.13) 


2.84 (0.14) 


0.65 


(0.15) 


1.35 (0.19) 


3C 305 


< 0.31 


3.82 (0.42) 


0.88 (0.09) 


1.89 (0.12) 


1.03 (0.39) 


2.93 (0.69) 


1.83 


(0.47) 


< 1.40 


3C459 


<0.25 


0.55 (0.03) 


<0.36 


0.37 (0.02) 


<0.29 


< 1.04 


< 


0.15 


<0.49 


4C 12.50 


2.31 (0.13) 


3.44 (0.39) 


0.96 (0.09) 


1.44 (0.09) 


<0.72 


< 1.41 


< 


0.93 


<2.63 


IC 5063 


4.63 (0.24) 


10.65 (0.44) 


< 1.8 


3.86 (0.49) 


< 3.12 


< 5.10 


< 


3.38 


<5.20 


OQ 208 


<0.38 


1.11 (0.09) 


<0.35 


0.49 (0.06) 


<0.75 


< 4.13 


< 


2.14 


< 2.76 


PKS 1549-79 


<0.46 


1.45 (0.12) 


0.63 (0.08) 


0.95 (0.09) 


< 1.29 


<2.23 


< 


1.32 


< 3.36 



Note. — H 2 mid-IR line fluxes (and lcr en'or in parenthesis) in units of 10 W m measured with Spitzer IRS. In case of 
non-detection, the 2a upper limit is indicated. The H 2 S(0) to S(3) lines are measured in the high-resolution modules (except 
for IC 5063, where the S(3) line was measured in the SL module because it falls out of the SH wavelength coverage), and the 
S(4) to S(7) lines in the low-resolution modules. Note that the S(5) and S(6) lines are severely blended with the PAH 6.2 fim and 
[ArII]A6.98 fim features. 



Table 5 







H 2 excitation diagram 


fitting for 


a thermal distribution: 


results 




Source 


Size" 


T(K) b 


o/p c 


W(H 2 ) d 


M(H 2 ) e 


L(H 2 ) f 




["] 


[K] 




[cm -2 ] 


[Mq] 


[L Q ] 






100 ( 0) 


1.587 


3.3E+21 


1.6E+09 (0.5) 


1.5E+07 


3C 236 


3.0 


241 (27) 


2.947 


1.9E+20 


9.5E+07 (4.0) 


8.2E+07 






1045 (80) 


3.000 


1.5E+18 


7.3E+05 (1.8) 


1.8E+08 






100 ( 0) 


1.587 


6.1E+21 


1.7E+09 (0.2) 


1.2E+07 


3C 293 


4.7 


285 (14) 


2.983 


1.9E+20 


5.3E+07 (1.1) 


7.4E+07 






1048 (47) 


3.000 


1.9E+18 


5.2E+05 (0.8) 


1.1E+08 






226(13) 


2.922 


<4.2E+20 


<1.0E+08 


5.3E+07 


3C 305 


4.7 


1253 (93) 


3.000 


1.3E+18 


3.1E+05 (0.5) 


1.2E+08 






100 ( 0) 


1.587 


<2.1E+22 


<8.3E+09 


1.1E+08 


3C 459 


1.4 


284 (32) 


2.983 


1.5E+20 


6.1E+07 (3.0) 


1.6E+08 






1500 ( 0) 


3.000 


1.6E+18 


6.3E+05 (1.1) 


7.9E+08 






100 ( 0) 


1.587 


2.2E+22 


3.8E+10 (0.2) 


3.7E+08 


4C 12.50 


4.7 


275 (16) 


2.978 


1.6E+20 


2.8E+08 (0.9) 


4.4E+08 






1500 ( 0) 


3.000 


5.5E+17 


9.2E+05 (1.5) 


8.3E+08 






100 ( 0) 


1.587 


3.9E+22 


7.3E+08 (0.9) 


4.7E+06 


IC 5063 


4.7 


210 (34) 


2.883 


1.4E+21 


2.6E+07 (1.9) 


8.8E+06 






1251 (103) 


3.000 


3.1E+18 


5.9E+04 (2.5) 


1.9E+07 






100 ( 0) 


1.587 


<9.2E+21 


<1.2E+09 


1.0E+07 


OQ208 


2.0 


186(29) 


2.784 


2.1E+21 


2.8E+08 (2.5) 


7.1E+07 






1500 ( 0) 


3.000 


2.5E+18 


3.4E+05 (0.4) 


2.6E+08 






100 ( 0) 


1.587 


<2.0E+22 


<7.7E+09 


8.4E+07 


PKS 1549-79 


1.9 


270(21) 


2.975 


5.1E+20 


2.0E+08 (0.7) 


3.3E+08 






1500 ( 0) 


3.000 


3.6E+18 


1.4E+06(0.2) 


1.4E+09 



a Assumed size (diameter) of the H 2 emitting source for the calculation of the column densities. 

b Fitted excitation temperature assuming LTE (see text for details). 

H 2 ortho-to-para ratio (fitted self-consistently to fulfill the LTE approximation). 

d H 2 model column densities. 

e Warm H 2 masses (and 1 <r uncertainties) derived from the fit. In case of non-detection of the S(0) line, a 2<r upper limit is 
quoted. 

f Total H 2 luminosity, summed over all rotational transitions (including those not observed with Spitzer, which includes a » 40% 
correction at 200 K). 



luminosity ratio that ranges from Lpah 11. 3/^24^1 = 4 x 10 2 
for 3C 293 to 8 x 10" 4 for 4C 12.50, which is similar to 
the ratios observed in the Ogle et al. (2010) sample of 3C 
radio-galaxies. This ratio is significantly lower than the 
SINGS star-forming galaxies, where the median ratio is w 0.1. 
3C459, PKS 1549-79 and 4C 12.50 are ULIRGs, but their 
ipAH 11. 3/^24 nm ratios are at the lower end of the observed val- 
ues (« 10~ 3 ). The mid-IR continuum of these radio-galaxies 
is therefore dominated by dust emission from the accretion 
disk or synchrotron emission from the jet, but not by star for- 
mation. The weakness of the PAH emission relative to the 
continuum emission is related to the weakness of the stellar 
UV radiation field (see Ogle et al. 2010, and sect. 6. 1 for a dis- 



cussion of this claim). Although PAH emission is expected to 
be diminished because of PAH destruction by hard UV or X- 
ray radiation close to the AGN Voit (1992), or because of gas 
suppression close to the nuclei, it is unlikely that their abun- 
dances are reduced significantly at galactic scales by AGN 
effects (Ogle et al. 2010). 

Silicate absorption at 9.7 and 18 /im is present in the low- 
resolution spectra of OQ 208 and 4C 12.50 (Willett et al. 
2010), as well as PKS 1549-79. OQ 208 is the only galaxy 
showing silicate emission (Willett et al. 2010). 



4.3. Ionic fine-structure lines 
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log 10 ( [Nelll]15.5fim / [Nell]12.8(jm ) log 10 ( [Nelll]15.5fim / [Nell]12.8(jm ) 

Figure 4. Mid-infrared line diagnostic diagrams for the logarithm of the model line flux ratios [NeV]/[NeII] (left) and [OlV]/[NeII] (right) against [NeIIl]/[NeII]. 
The observed values for the 8 Hi-outflow radio-galaxies are compared to shock and photoionization grids of models. The shock model grid is from the Allen et al. 
(2008) MAPPINGS III shock code for combined precursor and shock models. The gas metallicity is solar and the pre-shock gas density is «h = 10 cirT 3 . The 
orange lines indicate constant shock velocities, from 200 (light orange) to 1000 km s (red), and the blue lines constant magnetic field strengths, from 10~ 3 (light 
blue) to 100 /J.G (purple). The deeper the color, the larger is the parameter. The grid of photoionization NLR models was computed with the CLOUDY code 
(Ferland et al. 1998) with a standard power-law spectral index a = — 1 .5 and varying ionization parameter (green lines, Uo = 1 — 10~ 2 ) and density rm (pink lines, 
«H = 10 3 — 10 6 cm -3 ). On the right plot, since the shock and photoionization models partially overlap, the constant density lines are dashed for clarity. 



Table 6 

H2 fluxes fitted with shock models: results for «h = 10 3 cm" 



Source 


V s a 
[kms -1 ] 


Mass flow b 
[Mq yr 1 ] 


tcool(50K) c 
[yrs] 


M(Hi) d 
[Mq] 


3C 236 


3 

23 


8.0E+04 
2.4E+03 


1.5E+05 
1.8E+04 


1.2E+K) 
4.4E+07 


3C 293 


5 

40 


4.2E+04 
8.3E+02 


8.1E+04 
7.8E+03 


3.4E+09 
6.4E+06 


3C 305 


11 

40 


7.3E+03 
3.7E+02 


4.4E+04 
7.8E+03 


3.2E+08 
2.9E+06 


3C 459 


4 

40 


2.6E+05 
1.8E+03 


1.1E+05 
7.8E+03 


2.9E+K) 
1.4E+07 


4C 12.50 


4 

23 


1.1E+06 
8.6E+03 


1.1E+05 
1.8E+04 


1.2E+H 
1.5E+08 


IC 5063 


4 

40 


1.9E+04 
7.7E+01 


1.1E+05 
7.8E+03 


2.1E+09 
6.0E+05 


OQ 208 


5 

40 


5.4E+04 
3.6E+02 


8.1E+04 
7.8E+03 


4.3E+09 
2.8E+06 


PKS 1549-79 


15 
40 


2.3E+04 
1.9E+03 


2.9E+04 
7.8E+03 


6.7E+08 
1.5E+07 



a Shock velocity. 

b Mass of gas that is traversed by the shock per unit time required to match 
the observed H2 line fluxes. 

c Gas cooling time, from the peak post-shock temperature down to 100 K. 
d Warm tb masses. 



The fluxes of the most commonly detected mid-IR ion- 
ized gas lines are listed in Table 9. The [Nell] 12.8/im 
and [Nelll] 15.55/im are detected at high signal to noise ra- 
tio (> 6) in all 8 radio-galaxies. The [Olv]25.89/im and 
[Fell]25.99/im pair is also very common, with at least one 
line of the doublet detected in 7/8 galaxies. The high ioniza- 
tion (E ion = 97.12 eV) [Nev] 14.32^im line is detected in half 
of the sample. Although listed as non-detected by Farrah et al. 
(2007) and Willett et al. (2010), we securely detect this line in 
4C 12.50 at 3.2 a. 

Three objects show detections in both [Nev] 14.32 and 
24.32 /im lines. The flux ratios of these two lines are consis- 
tent with electron densities of n e m 3 x 10 3 cm -3 for IC 5063 



Table 7 

H2 fluxes fitted with shock models: results for «h = 10 4 cm 



Source 


V s a 
[kms -1 ] 


Mass flow b 
[M yr 1 ] 


t cool (50 K) c 
[yrs] 


M(H 2 ) d 
[Mq] 


3C 236 


4 
29 


2.2E+05 
2.8E+03 


2.5E+04 
1.8E+03 


5.6E+09 
5.0E+06 


3C 293 


8 

29 


4.0E+04 
1.5E+03 


1.1E+04 
1.8E+03 


4.4E+08 
2.7E+06 


3C 305 


5 

32 


5.5E+04 
1.3E+03 


1.8E+04 
1.5E+03 


1.0E+09 
1.9E+06 


3C 459 


8 

36 


8.3E+04 
6.2E+03 


1.1E+04 
1.2E+03 


9.2E+08 
7.5E+06 


4C 12.50 


4 
13 


1.8E+06 
3.2E+04 


2.5E+04 
5.3E+03 


4.4E+10 
1.7E+08 


IC 5063 


4 
38 


3.1E+04 
1.5E+02 


2.5E+04 
1.1E+03 


7.8E+08 
1.6E+05 


OQ 208 


3 

40 


1.6E+05 
1.6E+03 


3.7E+04 
1.0E+03 


5.7E+09 
1.7E+06 


PKS 1549-79 


8 

40 


1.7E+05 
1.1E+04 


1.1E+04 
1.0E+03 


1.8E+09 
1.2E+07 



a Shock velocity. 

b Mass of gas that is traversed by the shock per unit time required to match 
the observed H2 line fluxes. 

c Gas cooling time, from the peak post-shock temperature down to 100 K. 
° Warm tb masses. 



and PKS 1549-79 17 , and n e < 10 3 cm" 3 for 3C 305 (see e.g. 
Figure 3 of Alexander et al. 1999). These densities are be- 
low the critical densities of both lines (see also Dudik et al. 
2009). Using the [Sill] 18.71 and 33.48 line ratios, we 
find n e « 2 x 10 2 cm" 3 for 3C 305 and n e w 3 x 10 3 cm" 3 
for IC 5063. The Neon and Sulfur line ratios give consistent 
estimates of electron densities, although star-forming regions 
may contribute to the [Sill] line emission, whereas the [Nev] 
line is principally emitted by AGN-heated regions. 

Some more unusual lines, rarely seen in AGN, and that are 
not listed in Table 9, are detected in the high-quality spec- 

17 Holt et al. (2006) derive n e sa 430 era" 3 for the NLR gas using the 
[SII](6731/6717A) line ratio. Holt et al. (2011) derive n e = 3 X 10 3 cm- 3 
for 4C 12.50 using a more sophisticated technique. 
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Table 8 

Fluxes of the main PAH complexes from Hi-outflow radio galaxies 



Source 


PAH 6.2 fim 


PAH 7.7 /jm 


PAH 11. 3 pxa 


PAH 12.6/tm 


PAH 11 pm 


3C 236 


<2.5 


<4.5 


3.7(0.1) 


3.8 (1.2) 


3.2(1.3) 


3C 293 


11.1 (0.6) 


55.7 (2.5) 


21.5 (0.6) 


10.4 (0.4) 


25.6 (3.2) 


3C 305 


15.1 (3.4) 


63.1 (8.3) 


19.2(0.1) 


16.1 (0.2) 


11.2(2.6) 


3C 459 


7.2(1.3) 


22.5 ( 1.3) 


9.9 (0.2) 


9.7 (0.2) 


5.7 (0.3) 


4C 12.50 


9.8 (2.3) 


38.3(1.6) 


10.2 (0.2) 


18.1 (0.4) 


<9.5 


IC 5063 


<43 


189 (1.9) 


104.9 (2.6) 


230.6 (0.8) 


53.3 (6.7) 


OQ208 


40.8 (2.9) 


51.5 (1.4) 


55.3 (2.8) 


<61 


101 (8) 


PKS 1549-79 


25 (3.3) 


187 (5.3) 


37.4 (0.3) 


31.7 (0.7) 


4.8 (0.6) 



Note. — PAH fluxes (and lcr error in parentheses) in units of 10 17 W m 2 measured with 
Spitzer IRS. In case of non-detection, the 2<r upper limit is indicated. 



Table 9 

Fine-structure line fluxes in Hi-outflow radio-galaxies 



Object 


[Aril] 


[Arm] 


[Siv] 


[Nell] 


[Nem] 


[Sill] 


[Oiv] 


[Fell] 


[Sill] 


[Sin] 




6.98 Atm 


8.99 fim 


10.51 (im 


12.81 /xm 


15.55 pm 


18.71 /im 


25.89 fim 


25.99 pm 


33.48 /xm 


34.81 pm 


3C 236 


< 1.6 


< 1.9 


< 0.13 


0.92 (0.06) 


0.44 (0.03) 


<0.14 


<0.14 


< 0.13 


< 0.09 




3C 293 


2.6 (0.3) 


0.70 (0.12) 


< 0.18 


4.07 (0.29) 


1.27 (0.09) 


<0.42 


0.27 (0.02) 


0.57 (0.05) 


1.32(0.12) 


5.55 (0.48) 


3C 305 


2.8 (0.4) 


< 1.5 


3.21 (0.24) 


4.15 (0.29) 


8.62 (0.59) 


4.11 (0.33) 


13.07 (1.15) 


0.90 (0.07) 


5.14(0.46) 


11.24(0.97) 


3C 459 


1.0(0.1) 


< 0.15 


0.11 (0.01) 


3.11 (0.22) 


1.00(0.07) 


0.39 (0.03) 


0.46 (0.04) 


0.58 (0.05) 


a 


a 


4C 12.50 


<2.7 


2.4 (0.2) 


<0.7 


5.45 (0.39) 


2.99 (0.21) 


<0.83 


2.88 (0.25) 


<0.21 


a 


a 


IC 5063 


4.4 (0.2) 


9.5 (0.3) 


37.60 (3.05) 


22.83 (1.62) 


56.85 (3.92) 


23.88 (1.93) 


69.32 (6.10) 


3.97 (0.33) 


31.89 (2.68) 


28.3 (1.1) b 


OQ 208 


1.2(0.2) 


< 1.8 


< 1.6 


3.33 (0.24) 


2.60(0.18) 


< 1.12 


1.21 (0.11) 


0.45 (0.04) 


< 1.4 b 


a 


PKS 1549-79 


3.2 (0.3) 


3.3 (0.2) 


1.45 (0.12) 


2.20 (0.17) 


5.45 (0.39) 


<0.6 


3.11 (0.27) 


<0.8 


a 


a 



Note. — Ionic line fluxes (and lcr error in parentheses) in units of 10~ 17 W irT 2 measured with Spitzer IRS. In case of non-detections, a 2cr upper limit is 
quoted. The [Aril] and [Arm] lines are measured in the low-resolution modules. All the other lines are measured in the high-resolution modules, except some of 
the [Sill] and [Sill] lines marked with a (b). 
11 Observed wavelength not visible in the IRS range. 
b Measured on the low-resolution LL module. 



Table 10 

High-ionization forbidden emission lines 



Object 


[NeVI]A7.65 pm 


[NeV]A14.32/jm 


[NeV]A24.32/tm 


3C 236 


< 1.6 


< 0.06 


<0.06 


3C 293 


<0.53 


<0.21 


<0.17 


3C 305 


0.92 (0.32) 


1.42(0.12) 


1.87 (0.15) 


3C 459 


<0.29 


<0.16 


<0.14 


4C 12.50 


2.66 (0.18) 


1.68(0.22) 


< 1.1 


IC 5063 


3.51 (0.12) 


22.70(1.86) 


17.41 (1.43) 


OQ 208 


< 1.2 


<0.8 


<0.2 


PKS 1549-79 


3.29 (0.21) 


3.35 (0.27) 


2.21 (0.18) 



Note. — Line fluxes (and lcr error in parentheses) in units of 10 W m 
measured with Spitzer IRS. In case of non-detection, a 2cr upper limit is 
quoted. The [NeVI]A7.65 pm line is measured in the SL module. The other 
lines are measured in the high-resolution modules. 

trum of 3C 305. The low-ionization (E[ on = 7.9 eV) [Fell] 
17.94, 24.52, and 25.99 fim lines are present, as well as the 
[Fem] 22.93/mi (£ ion = 16.19 eV). The [Arv] 13. 10/im line is 
also weakly detected (see the SH spectrum on Figure 9). 

In such objects, the excitation of the fine-structure 
lines can arise from photo-ionization by the AGN (e.g. 
Veilleux & Osterbrock 1987; Meijerink et al. 2007) or by 
shocks (e.g. Dopita & Sutherland 1996). In practice, dis- 
tinguishing between the two excitation mechanisms is diffi- 
cult, since the line flux ratios can be differentially affected by 
the structure (homogeneous or clumpy) of the emitting gas, 
and since the model predictions for fine-structure line flux ra- 
tios are degenerate for both excitation mechanisms (see also 
Groves et al. 2006). 

In Sect. 6.1, we will show that the warm IT emission is 



likely powered by shocks. Is it also the case for the ob- 
served ionized line emission? In Figure 4, we compare the 
observed values of three fine-structure line flux ratios to the 
predictions of shock and photo-ionization models. The grid 
of fast (> 200 km s" 1 ) radiative shocks is from Allen et al. 
(2008). The pre-shock gas density is «h = 10 cm" 3 . Note 
that augmenting the density by a factor 10 reduces these line 
flux ratios by a factor 2 - 3 on average, thus shifting the mod- 
els to the bottom left of Figure 4 by 0.3 - 0.5 dex. The 
photo-ionization models are computed with the CLOUDY 
code (Ferland et al. 1998), with a setup updated from the nar- 
row line region (NLR) models presented in Ferland & Netzer 
(1983). We used a dusty gas slab model at constant pressure, a 
power-law ionizing continuum with f v oc v and no cut-off 
or extinction. We varied the hydrogen density of the gas from 
10 3 to 10 6 cm -3 and the ionization parameter 18 from Uo = 10~ 2 
to 1 . Note that these models are simplistic: they do not include 
geometrical and inhomogeneity effects, as well as metallicity 
differences (assumed to be solar here). A detailed model of 
the mid-IR emission of these sources is beyond the scope of 
this paper. 

The [Nem] 15.55/xm / [Nell] 12.8/xm line flux ratios are in 
the range of 0.3-2.4, much higher than those observed in 
pure starburst galaxies (0.05-0.2, e.g. Bernard-Salas et al. 
2009). These ratios are consistent with shock models for v < 
250 km s" 1 and «h < 10 2 cm -3 , the higher values of [Nem] / 
[Nell] m 2 requiring very low gas densities («h < 1 cm" 3 ). 

18 Uo = Si,/(n e c) is a direct measure of the ionization state of the gas, 
where SV is the flux of ionizing photons, n e the electron density at the inner 
face of the slab of gas, and c the speed of light. 
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The range of observed [Nelll]/[Nell] ratios is also consis- 
tent with photo-ionization models. This ratio is weakly de- 
pendent on the ionization state of the gas, and sensitive to 
the gas density. The observed [Nelll]/[Nell] ratios require 
Kb > 10 3 cm" 3 . 

[Nev] 14.32 /im being a high ionization potential line 
(£ ion = 97.1 eV), the [Nev] 14.32/[Nell] 12.8 line ratio is 
very sensitive to the ionization conditions, and therefore a 
strong indicator of the AGN contribution (e.g. Sturm et al. 
2002). The shock radiative precursor, which ionizes the 
gas ahead of the shock front by UV and soft X-ray pho- 
tons, adds an emission component of a highly ionized gas. 
However, for the sources containing a powerful AGN, the 
high values of the observed [Nev]/[Nell] ratios (left panel 
in Figure 4) rule out shocks as the dominant excitation 
mechanism to power the [Nev] luminosity. The [Nev] 
line primarily arises from photo-ionization of the gas by 
the AGN. According to this grid of CLOUDY models, the 
observed [Nev] 14.32/[Nell] 12.8 ratios are consistent with 
n H > 10 3 cnT 3 andt/o <0.1. 

The right panel in Figure 4 uses the lower ionization (£ ion = 
54.93 eV) [Orv]25.89/im line. In this case, the shock and 
photo-ionization model grids cover similar regions of the dia- 
gram, making it difficult to distinguish between the two exci- 
tation mechanisms. The observed mid-IR line ratios are both 
consistent with high- velocity shocks (v > 200 km s" 1 ) and 
photo-ionization, with a preference for shocks at low values 
of [Nem]/[Neii]. Some sources (3C 459, 3C 293 and 3C 236) 
require higher pre-shock densities («h = 10 2 - 10 3 cm -3 ) to be 
consistent with shock models. 

5. KINEMATICS OF THE MULTIPHASE GAS 

In this section we describe and compare the kinematics of 
the mid-IR H2 and ionized gas lines (from the IRS data and 
optical data from the literature) to the kinematics of the out- 
flowing Hi gas. 

5.1. Analysis of the line profiles 

We fitted all the mid-IR lines with Gaussians on the dust- 
free high-resolution spectra (sect. 3.2). This allows us to care- 
fully extract line profiles, especially when those are located 
close to underlying PAH structure, as it is the case for the 
[Nell] 12.8/im line for instance. We pay a special attention 
to the fitting of the [Nell] by checking that the PAH 12.6 /mi 
has been correctly removed by the PAHFIT modelling, since 
the PAH features at 12.62 and 12.69 /im (Hony et al. 2001) 
may contribute to an apparent blue-shifted wing of the line. 
We validated our subtraction method by simulating the pres- 
ence of an artificial Lorentzian complex in the 12.6-12.9 /im 
range of intensity compatible with the observed PAH 6.2 and 
7.7 /im features. We find that the strength of the simulated 
12.6 /im feature measured with PAHFIT is at least a factor of 
5 smaller than the strength of the [Nell] blue wing. In ad- 
dition, the strength of the 12.7 /im PAH complex correlates 
with the 6.2 /im feature, which is another indication that our 
measurement of the 12.7 /im PAH complex is correct. 

Figures 10 and 11 show a detailed view of the individual 
lines. The blended lines, like [Oiv]A25.89 and [Fell] A 25. 99, 
and the lines exhibiting wings (see sect. 5.4) were fitted with 
a sum of two Gaussian components. 

We assume a Gaussian decomposition of the instrumental 
profile (shown on each of the lines in Figures 10 and 11), and 
we derive the intrinsic linewidth (FWHM w,- in km s -1 ) from 



Table 11 

Intrinsic H2 linewidths in Hi-outflow radio-galaxies 





H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 




28.22 fim 


17.03 /im 


12.28 


9.66 


3C 236 


<50() 


548 (59) 


< 500 


< 500 


3C 293 


743 (79) 


485 (54) 


524 (63) 


586 (68) 


3C 305 




482 (54) 


525 (63) 


566 (66) 


3C 459 




432 (50) a 




751 (81) 


4C 12.50 


906 (91) 


< 500" 


529 (63) 


< 500 


IC 5063 


709 (76) 


721 (71) 


< 500 




OQ 208 




< 500 






PKS 1549-79 




585 (62) a 


< 500 


< 500 


Note. — H 2 


FWHM in units of km s 1 measured by single Gaussian fit 



ting. A Gaussian decomposition from the instrument profile is assumed to 
derive the intrinsic linewidths (\Vj, see sect. 5.1 for details). 
a 3C 459, 4C 12.50 and PKS 1549-79 possibly exhibit a blue wing in the 
H 2 S(l) line, which is not included in the quoted linewidth here because this 
detection is marginal. 

w; = (Wm ~~ vv irs) 1/ ' 2 ' wner e w m is the FWHM in km s" 1 di- 
rectly measured on the observed spectrum, and wirs = c/TZ is 
the velocity resolution of the IRS. c is the speed of light, and 
1Z = 600 ± 72 is the spectral resolution of the high-resolution 
module, corresponding to wirs = 500 ± 60 km s" 1 . Measure- 
ments of calibration targets indicate that the resolution of the 
high-resolution module of the IRS is constant over its wave- 
length coverage (Dasyra et al. 2008). The error on w m is com- 
puted as e„, = (e/ + e 2 Rs)'^ 2 > where eiRs = 60 km s" 1 is the error 
on the instrumental resolution and ej the error of the mea- 
surement, estimated from two Gaussian fits of the line profile, 
one with the upper flux values (adding flux uncertainties), one 
with the lower values (subtracting uncertainties). A line is 
considered resolved when w m - e m > wirs + £irs = 560 km s -1 . 
Unresolved lines are assigned a conservative upper limit of 
500 km s" 1 on their intrinsic linewidths. 

5.2. Warm Hi gas 

The intrinsic linewidths (w,) of the H2 rotational lines are 
listed in Table 1 1 . According to the criterion defined above, 
more than 60% of the detected H2 lines are resolved by the 
Spitzer IRS, with FWHM up to rj 900 km s" 1 . All the sources, 
except OQ 208, have at least one broad H2 line with intrinsic 
FWHM> 430 km s" 1 . 3C 293, 3C 305, 3C 459, 4C 12.50, and 
IC 5063 show two or more resolved H2 lines. These results 
show that the warm H2 gas is very turbulent in most of the 
sources. As we will show in sect. 6.2, a fraction of the warm 
H2 emission is at velocities larger than the escape velocity. 

Three of the sources, 3C 459, 4C 12.50 and PKS 1549-79 
exhibit asymmetric H2 line profiles, with blue-shifted wings 
(Figures 6, 10 and 1 1). However, the detection of these wings 
is tentative, at a 2.5- 2.8 a significance for the S(l) line, and 
even weaker for the S(2) and S(3) lines. Therefore we do not 
attempt to fit these wings, and we do not include them in the 
FWHM quoted in Table 1 1 . Deeper observations at higher 
spectral resolution are needed to confirm the blue-shifted H2 
emission and measure what fraction of the the warm H2 gas 
is participating to the outflow, and to derive the outflow pa- 
rameters (warm H2 mass in the outflowing component, mass 
outflow rate, etc.). 

5.3. Ionized gas 

Table 12 lists the intrinsic linewidths of the fine-structure 
lines. Remarkably, the [Nell] 12.8/im and [Nelll]15.5/im 
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Table 12 

Intrinsic linewidths of ionic fine-structure lines in Hi-outflow radio-galaxies 



Object 


[SIV] 


[Nell] 


[Nev] 


[Nein] 


[Sill] 


[Nev] 


[Oiv] 


[Fell] 


[Sill] 


[Sin] 




10.51 /im 


12.81 Aim 


14.32 (Um 


15.55 fj,m 


18.71 /im 


24.32 Atm 


25.89Atm 


25.99 Atm 


33.48 Atm 


34.81 A*m 


3C 236 




1239(109) 




875 (86) 














3C 293 




861 (84) 




877 (86) 






723 (83) 


595 (70) 




779 (85) 


3C 305 


709 (79) 


864 (84) 


611 (72) 


817 (82) 


543 (64) 


<500 


488 (62) 


< 500 


565 (67) 


746 (82) 


3C 459 


573 (68) 


775 (77) 




<500 


794 (125) 




< 500 


789 (86) 






4C 12.50 




1033 (95) 


852 (92) 


838 (83) 






767 (87) 


< 500 






IC 5063 


< 500 


408 (49) 


< 500 


<500 


632 (71) 


<500 


< 500 


661 (75) 


<500 




OQ 208 




< 500 




895 (87) 






781 (89) 


< 500 






PKS 1549-79 


1355 (132) 


848 (83) 


967 (89) 


1233 (105) 




834 (89) 


857 (95) 









Note. — FWHM in units of km s 1 measured by single Gaussian fitting. A Gaussian decomposition from the instrument profile is assumed to 
derive the intrinsic linewidths (wj, see sect. 5.1 for details). 



lines, the brightest among detected fine-structure lines, are 
spectrally resolved in 6/8 objects, with intrinsic FWHM up 
to w 1250 km s" 1 (velocity dispersions up to rj 540 km s" 1 ). 

In five sources (3C 236, 3C 293, 3C 459, 4C 12.50 and 
PKS 1549) the [Nell], and to a lesser extent [Nelll] and 
[Nev], line profiles are asymmetric. We securely detect 
blue-shifted (and in some cases red-shifted) wings, up to 
3000 km s" 1 wide, underlying a stronger and narrower peak 
component centred very close to the optical systemic veloc- 
ity. As explained in sect. 5.1, we carefully checked that the 
[Nell] 12.8/^m line wings cannot be ascribed to an underlying 
PAH spectral feature. This result is confirmed by the fact that 
in most cases the [Nelll]15.5/im line profile, though noisier 
but located in a PAH-free part of the spectrum, is consistent 
with the [Nell]12.8/im line profile. The kinematic proper- 
ties of the broad [Nell] 12.8/.tm features are listed in Table 13. 
Such wings in the Neon lines have been reported in 4C 12.50 
(and in other ULIRGs) by Spoon & Holt (2009). 

To characterize the asymmetry and study the kinematics of 
the outflowing ionized gas, we decomposed the asymmetric 
line profiles into two Gaussian components, which provides 
a satisfactory fit. The resulting fit and the broad component 
are shown on Figures 10 and 11. The detailed results of 
the fits are listed in Table 13. The intrinsic FWHM of the 
broad [Nell]12.8Avmi components range from « 730 km s" 1 
(PKS 1549-79) to « 2300 km s" 1 (3C 293). In three sources 
(3C 459, 4C 12.50, and PKS 1549-79), the broad component 
is clearly blue-shifted with respect to the systemic velocity, 
with velocity shifts up to w 800 km s" 1 . For 3C 236, the en- 
tire [Nell] 12.8/im line and the broad Hi absorption profile are 
redshifted by ps 500 km s -1 with respect to systemic, whereas 
the H2 line is centred on the systemic velocity. This is a clear 
case where a large fraction of the ionized and Hi gas is in an 
outflow, but the H2 gas is not. 

5.4. Hi gas 

In this section we use the WSRT Hi absorption profiles 
(Morganti et al. 2005b) to quantify the kinematics and ener- 
getics of the Hi gas. We recall that this Hi absorption, detected 
against the radio continuum, is spatially unresolved, except 
for IC 5063 and 3C 305. The profiles all exhibit a deep and 
narrow (FWHM<200 km s" 1 ) component that traces quies- 
cent gas (likely located in a large scale disk), plus shallow and 
broad components (with FWHM up to 970 km s" 1 ), mostly 
blueshifted, that indicate outflowing gas (e.g. Morganti et al. 
2003). Two or three Gaussians were necessary to accurately 
fit the profiles (see Figure 15). The Gaussian fitting parame- 
ters and energetics of the Hi gas are summarized in Table 14. 



At each given velocity, the optical depth (r) is calculated 
from e~ T = 1 -S a bs/S C ont , where S a bs is the absorption flux 
and Seem is the underlying radio continuum flux. Then, the Hi 
column density is derived from: 



Nhi 



= 1.822 x 



1 spin 

10 3 K 



r(v)dv , 



(2) 



where v is the velocity and r sp ; n is the spin temperature of the 
Hi gas, assumed to be 1000 K, given that the Hi gas is heated 
by X-ray photons and shocks (see e.g. Bahcall & Ekers 1969; 
Maloney et al. 1996). Note that this value of r sp i n is an order 
of magnitude, and could well be within a 10 2 - 10 4 K range. 
The mass of Hi gas is derived from Nhi and the size of the Hi 
absorption region (see sect. 2 and Table 2). 

The turbulent kinetic energy of the Hi gas (associated with 
the velocity dispersion of the gas) is estimated from E^ b ' = 
3/2Mhict 2 ii , where ohi = FWHM/2.36 is the velocity disper- 



7 (turb) 



takes into account 



sion of the Hi gas. The factor of 3 in 
the three dimensions. 

The bulk (radial) mechanical energy of the entrained gas 
in the outflow is calculated directly from the broad Hi profile 
(after subtraction of the narrow component) by integrating the 
following quantity relative to the systemic velocity (v sys ): 



rbulk 

=392_ = 1.23E4 



In- 



Set 



-V 



dv 



lkpcl0 3 K7 5 con t-5abs V10 3 km/s/ lkm/s 

(3) 

where r is the Hi outflow radius (see Table 2 and 
Morganti et al. 2005b). 

For comparison, one can assume an outflow with a constant 
velocity and mass flow M, and derive the energy loss rate as 
(Heckman 2002): 



£ = 8.1 x 10 8 C 



Nhi 



4tt 1 kpc 10 21 cm" 2 \350kms-' 



(4) 



where C/ is the covering fraction and ft the opening solid an- 
gle in which the gas the outflowing from a radius r. 

The energetic quantities defined above are listed in Ta- 
ble 14. The conversion from the kinetic energies to the kinetic 
luminosities are done via a timescale defined with the radius 
r and the outflow velocity (or velocity dispersion). The bulk 
(radial) Hi kinetic luminosities and energy loss rates are larger 
than the turbulent (velocity dispersion only) kinetic luminosi- 
ties, except for OQ 208. The size of the OQ 208 radio source 
is small and the Hi absorption spectrum exhibits three broad 
components, which makes the identification of the outflow 
component difficult. The bulk Hi kinetic luminosities repre- 



Spitzer IRS spectroscopy of Hi-outflow radio galaxies 



13 



Table 13 

Ionized gas kinematics: fit results for mid-IR lines with broad wings 



Line 


Comp. 






obs. FWHM w m 


int. FWHM Wj 


Flux 






[/an] 


[kms 4 ] 


[kms 4 ] 


[kms 4 ] 


[10-"Wm- 2 ] 










3C 236 






[Nell] 


n 


12.792 


-507 ± 33 


981 ±96 


Oil 1 O O 

844 ± 82 


0.55 ±0.04 


[Nell] 


b 


12.796 


AO A 1 ^ O ^ 

—424 ± 121 


21 11 ± 201 


21 14 ± 195 


0.38 ± 0.03 










3C 293 






[Nell] 


n 


12.818 


98 ±26 


779 ± 85 


597 ±65 


2.66±0.19 


[Nell] 


b 


12.835 


494 ±51 


2338 ±215 


2284 ± 209 


1.91 ±0.16 










3C 459 






[Nell] 


n 


12.824 


227 ±21 


804 ± 86 


630 ± 67 


2.50±0.18 


[Nell] 


b 


12.797 


-372 ± 89 


1721 ±163 


1647 ±155 


0.80 ±0.06 










AC 12.50 






[Nell] 


n 


12.810 


-88 ±15 


786 ±85 


606 ± 65 


3.35 ±0.24 


[Nell] 


b 


12.779 


-812 ±99 


2097 ± 194 


2037 ± 189 


2.46 ±0.20 










PKS 1549-79 






[Nell] 


n 


12.814 


-5 ±12 


793 ± 83 


616 ±66 


1.86±0.13 


[Nell] 


b 


12.775 


-906 ±131 


885 ± 142 


730 ±82 


0.48 ± 0.04 



Note. — Results from a two-component Gaussian decomposition of the [NeII]12.8/tm emission 
line profiles into narrow (n) and broad (b) components. Ac is the central wavelength of the com- 
ponent and corresponds to a velocity shift v c of the component with respect to systemic velocity. 
The statistical error on the velocity shift does not include corrections for the uncertainty in the 
wavelength calibration, which is 1/5 of a resolution element (100 km s~'). w m (respectively w,) 
is the observed (resp. intrinsic) FWHM before (resp. after) removal of the instrumental profile in 
quadrature (see sect. 5.1 for details). 



Table 14 

Hi gas kinematics and energetics: broad Hi component associated with the neutral outflow 



Object 


FWHM a 
[km s _1 ] 


Velocity b 
[kms" 1 ] 


timescale c 
[yrs] 


r d 


w (tmb) e 
c kin 

[Joules] 


„(bulk) f 

c kin 

[Joules] 


T (lurb) g 

[L Q ] 


j {bulk) h 
[L ] 


£(HI) ' 
[L ] 


3C 236 


932 


29367 


6.5E05 


0.0033 


5.8E49 


2.5E50 


7.2E09 


1.6E10 


2.0E10 


3C 293 


852 


13156 


2.0E06 


0.0038 


3.9E49 


3.9E50 


1.6E09 


1.1E10 


1.4E10 


3C 305 


689 


12402 


3.9E06 


0.0023 


1.2E49 


3.9E50 


2.5E08 


1.2E10 


5.9E08 


3C 459 


411 


65934 


3.3E06 


0.0005 


1.6E48 


1.2E50 


4.1E07 


4.0E09 


3.8E08 


4C 12.50 


970 


36617 


3.3E05 


0.0017 


5.1E48 


2.7E49 


1.3E09 


5.7E09 


2.1E09 


IC 5063 


133 


2785 


1.1E06 


0.0120 


2.3E48 


1.4E50 


1.7E08 


8.2E09 


3.2E09 


OQ 208 


482 


22568 


1.6E04 


0.0057 


5.0E48 


2.0E47 


2.5E10 


6.5E08 


3.4E08 


PKS 1549-79 


329 


45553 


3.9E06 


0.0200 


5.9E48 


3.3E51 


1.2E08 


1.0E11 


1.2E10 



a Full width at half maximum of the broad Hi component (fitted with a Gaussian). 
b Central velocity of the broad Hi component. 

c Outflow timescale derived from the radius r and the outflow velocity. 
d Optical depth of the broad Hi component for each velocity components. 
e Turbulent kinetic energy computed from 3/2M(HI)<r 2 . 

f Bulk kinetic energy computed by integrating the Hi profile (see sect. 5.4 for details). 
g Turbulent kinetic luminosity. 



Bulk Hi kinetic luminosity computed from Eq. 3 
1 Hi kinetic energy loss rate computed from Eq. 4 

sent about a few percent of the Eddington luminosities. This 
is two orders of magnitude larger than the ionized outflow ki- 
netic power (Holt et al. 2006; Morganti et al. 2007, 2010). We 
will further discuss these results in sect. 6.2 when examining 
the transfer of kinetic energy from the jet to the Hi and H2 
gas. 

5.5. Comparison of the multiphase gas kinematics 

In Figure 5, we compare the measured FWHM of the H2 
S(l) and [Nell]A12.8/im lines. Except for IC 5063 and 
OQ 208, the [Nell] is systematically broader, by a factor of 
40%, than the 17 /im H2 S(l) line. This is also the case for the 
other ionic fine-structure lines: the velocity dispersion of the 
ionized gas is generally larger than that of the warm Ho gas. 

Figure 6 compares the kinematics of the ionized, atomic 
and warm H 2 gas. In 3C 236, 3C 293, 3C 305, 3C 459, 
4C 12.50 and PKS 1549-79, the blue-shifted, broad Hi ab- 



sorption profile matches the blue wing of the [Nell] A 12.8/mi 
line. The ionized and atomic gas are well coupled dy- 
namically and outflowing at comparable velocities. Most 
of these galaxies are also known to exhibit outflows from 
optical spectroscopy: blue-shifted 5007 A[Olll] line emis- 
sion have been reported in 3C 293 (Emonts et al. 2005), 
3C 305 (Morganti et al. 2005b), IC 5063 (Morganti et al. 

2007) , PKS 1549-79 (Tadhunter et al. 2001; Holt et al. 2008), 
4C 12.50 (Holtetal. 2003, 2008) and 3C 459 (Holtetal. 

2008) . For these three targets, the [Olll] velocity widths and 
shifts with respect to the systemic velocity are very close to 
our measurements of the [Nell] A 12.8/im line. Interestingly, 
in PKS 1549-79, the [Nev]A 14.3/xm line is blueshifted by 
500-700 km s~\ at the position of the [Nell] and [Nelll] 
blue wings. This increase of blueshift with increasing ioniza- 
tion has been observed in ULIRGs (Spoon & Holt 2009), and 
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Figure 5. [Nell] A 12.8/mi versus H 2 S(l) 17/im linewidths (FWHM in 
km s _I ) measured from the Spitzer IRS spectra with the high-resolution mod- 
ule (w m in sect. 5.1). The linewidths are intrinsic, i.e. they have been cor- 
rected for the instrumental Gaussian profile. 

suggests that the outflow speed decreases with distance to the 
central ionizing source. 

On the other hand, in general, the kinematics of the warm 
H2 gas does not follow that of the ionized or atomic gas. In 
the cases where a blue wing is tentatively detected on the H2 
S(l) line (for 4C 12.50, 3C 459 and PKS 1549-79), the ve- 
locity extent of this H2 wing is smaller than that of the [Nell] 
line. If some of the H2 gas is entrained in the flow, its ve- 
locity is at least a factor of 2 to 3 smaller than that of the 
ionized and atomic gas, based on the comparison between the 
velocity extent of the blue-shifted H2 S(l) signal, and the one 
observed for the optical and mid-IR ionized gas lines. Note 
that the broad molecular line - e.g. CO(1-0)- components 
detected in composite sources like Mrk 231 are at a 1 -5% 
level of the peak line flux. The limited sensitivity and spectral 
resolution of the Spitzer IRS do not allow us to securely de- 
tect such broad wings in the H2 lines at comparable levels. In 
addition, if the molecular gas in these galaxies lies in a very 
turbulent rotating disk (as confirmed by our VLT/SINFONI 
observations of 3C 326, Nesvadba et al. 201 la), H2 line wings 
in spectra integrated over the whole galaxy could be difficult 
to to identify on top of the fast rotation-velocity field. 

6. DISCUSSION 

6. 1. Which dominant mechanism powers the H2 emission ? 

Following Ogle et al. (2010), we use the H2 to PAH lumi- 
nosity ratio, plotted on Figure 7, to investigate the contribu- 
tion of UV photons to the total heating of the H2 gas. In 
"normal" star-forming galaxies and dwarfs, the tight corre- 
lation between the H2 and PAH luminosities, and the cor- 
respondence of the H2-to-PAH luminosity ratio with PDR 
models, suggest that most of the H2 line emission is pow- 
ered by stellar UV photons (e.g. Rigopoulou et al. 2002; 
Higdon et al. 2006; Roussel et al. 2007). Strikingly, all of 
the Hi-outflow radio-galaxies have a L(H2)/L(PAH7.7) ratio 
larger (by at least a factor of 4) than the median value ob- 
served for star-forming galaxies (0.0086 for the SINGS galax- 
ies, Roussel et al. 2007). These H2-bright galaxies fall in the 
MOHEGs (molecular hydrogen emission galaxies) category 
(Ogle et al. 2010), defined as L(H 2 S(0)-S(3))/L(PAH7.7) > 
0.04. Using the Meudon PDR code (Le Petit et al. 2006), 



we computed the H 2 S(0)-S(3) to PAH 7.7 pm flux ratio as 
a function of the Guv/"h, where Guv is the intensity of 
the UV radiation field (in Habing units) illuminating a slab 
of gas of hydrogen density «h- The H2-to-PAH ratio is 
high for small values of Guv/«h- Exploring densities from 
«h = 10 2 to 10 4 cm" 3 and Guv = 1 — 10 4 , we find a maxi- 
mum H2-to-PAH flux ratio of 4 x 10~ 2 , which is in agree- 
ment with the Kaufman et al. (2006) models. This value is 
precisely the limit chosen empirically by Ogle et al. (2010) to 
define MOHEGs. All of the Hi-outflow radio-galaxies have 
L(H 2 )/L(PAH7.7) > 0.04 (see Figure 7), showing that their 
H2 emission cannot be accounted only by UV heating. 

Interestingly enough, Figure 7 shows that the sources hav- 
ing the largest H2 to PAH luminosity ratio have moderate IR 
luminosities (L^A/xm ~ 10 9 L Q ) and are jet-dominated, suggest- 
ing that the jet plays an important role in powering the H2 
emission. At larger IR luminosities, the sources become dom- 
inated by star formation, hence the more moderate H2 to PAH 
ratios. We do not find any extreme H2 to PAH luminosity ra- 
tios (> 1) at large IR luminosities (£24^,,, ^ 1O U L0). This 
may be due to the fact that all of these sources are at low (< 
0.3) redshifts, so the most luminous quasars are not included 
in this sample. We note that the predominately radio-quiet 
Higdon et al. (2006) sample of ULIRGs (0.02 < z < 0.93) and 
the Schweitzer et al. (2006) sample of QSOs (z < 0.3) have H 2 
to PAH luminosity ratios significantly higher than the SINGs 
star-forming galaxies. Although their H2 to PAH luminos- 
ity ratios are compatible with the range of ratios predicted by 
PDR models, other studies (based on extinction or H2-to-CO 
line ratios), suggest that additional sources of H2 excitation 
(other than UV photon heating) are required to explain the 
H2 emission (e.g. shocks, see Zakamska 2010; Fioletetal. 
2010). 

AGN X-ray heating is potentially a source of H2 excita- 
tion (e.g. Maloney et al. 1996). Assuming a characteristic H2 
gas temperature of 200 K and that all the X-ray flux from the 
AGN is absorbed by the molecular gas, Ogle et al. (2010) es- 
timated that the maximum H2-to-X-ray luminosity ratio is of 
the order of L(H 2 0-0S(0)-S(3))/L x (2- lOkeV) < 0.01. We 
use the c08.01 version of the CLOUDY code, last described 
by Ferland et al. (1998), to compute this ratio for a cloud of 
density «h = 10 4 cm" 3 and a total ionizing photon luminosity 
corresponding to L x (2- lOkeV) = 10 43 erg s , which is the 
average X-ray luminosity measured for our sample. The cov- 
ering and filling factors are set to unity. For this set of parame- 
ters we findL(H 2 0-0S(0)-S(3))/L*(2- lOkeV) = 7 x 10" 3 , 
which is in reasonable agreement with the upper limit of 0.01 
given above. All Hi-outflow radio-galaxies that have X-ray 
measurements have an observed H2-to-X-ray flux ratio above 
that limit (see last column of Table 1), which shows that the 
X-ray heating cannot be the dominant powering source of the 
H2 emission in Hi-outflow radio MOHEGs, except perhaps in 
IC 5063, where the H^-to-X-ray luminosity ratio is close to 
0.01. 

Another source of heating for the warm H2 gas is cosmic 
rays, as proposed by Ferland et al. (2008) to explain the mid- 
IR H2 emission observed from the Ha-emitting filament in 
the Perseus A cool-core cluster (Johnstone et al. 2007). Fol- 
lowing, the discussion by Ogle et al. (2010), we find the ob- 
served H2 S(0)-S(3) line cooling (i.e. the H2 luminosity per 
H2 molecule, using the modelled H2 masses given in Ta- 
ble 6) could be balanced by cosmic ray heating if the cosmic 
ray ionization rate (per H) is ( H = 2 x 10~ 14 -3 x 10~ 13 s -1 , 
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Figure 6. Comparison of the normalized Hi absorption (data from Morganti et al. 2005b) (shown on a positive scale) and the Spitzer IRS H2 S(l), [Nell] A 12.8 fim 
and [NeV]A 14.3 fim (when detected) emission profiles. The dotted black line is the shallower broad Hi component associated with the outflowing gas, after 
removal of the deep, narrower absorption peak. The optical systemic velocity is indicated by a vertical red dashed line. For clarity, the Hi spectra of 3C 236, 
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Figure 7. Ratio of the mid-IR H2 line luminosities (summed over S(0) to S(3)) to the PAH U/im emission vs. narrow-band 24/^m continuum luminosity. This 
ratio indicates the relative contribution of mechanical heating (shocks) and star-formation (SF) power. The red pentagons are the nearby radio galaxies with fast 
(> 1000 km/s) Hi outflows. The orange triangles are the H2-bright radio-galaxies presented by Ogle et al. (2010), augmented by the compact symmetric objects 
observed by Willett et al. (2010). The green ellipses are elliptical galaxies from the Kaneda et al. (2008) sample. The purple triangles are the H2-bright cool-core 
galaxy clusters from Donahue et al. (2011). These H2-luminous galaxies stand out above SF and AGN galaxies from the SINGS survey (open diamonds and 
blue circles; data from Roussel et al. 2007). The black dotted line shows the median value of the H2-to-PAH luminosity ratio for the SINGs SF galaxies. The 
H2 emission in these sources cannot be accounted for by UV or X-ray photon heating. The blue dashed line shows the upper limit given by the Meudon PDR 
models (Le Petit et al. 2006) obtained for nu = 10 4 cm -3 and Guv = 10. For comparison, a few other types of H2-luminous galaxies are shown (black crosses): 
the Stephan's Quintet (SQ) and Taffy galaxy collisions (data from Cluver et al. 2010, and Peterson et al., sub.), ZW 3146 (Egami et al. 2006) and Perseus A 
(Johnstone et al. 2007) clusters, and the NGC 6240 merger (Armus et al. 2006). The black squares indicate averaged values for the Higdon et al. (2006) sample 
of ULIRGs. The black star is an averaged value of the Schweitzer et al. (2006)sample of QSOs. 



Spitzer IRS spectroscopy of Hi-outflow radio galaxies 



17 



"en 





''' ^ 4C12.50 /" 


jk 3C459 x PKS1 549 / 

/ 

/ 

/ 


/ _ 

/' 

/ / 
* 3C236 




/ m 3C293 / „ OQ208 


; / *3C305 , 






G)/ 


~ if' 




/// /// 




' i>' 


■ \ ' 1 



10 7 10 s 10 9 10 10 10 1 ' 

HI turbulent kinetic luminosity [L ] 





10° 




' / 1 ' ' ' * 401'2.50 
* 30459 


.,■ PKS1 549" 


_] 

^ 

"cn 

O 
c 

E 












4 3C236 






10 8 




» 0Q208 / j/3C293 










,' m 3C305 




(/) 










o" 






A/' 




t/T* 




"? / 






X 


10 7 


// >' 


// / 

&f M IC5063 

1 ■ ... 


. , , . i 




10 s 


10 9 10 10 


10" 








HI bulk kinetic luminosity [L ] 





Figure 8. Top: H2 luminosity (summed over the S(0)-S(3) lines) against 
the Hi turbulent kinetic luminosity, derived from the linewidth of the broad 
Hi absorption component, associated with the outflowing gas. Bottom: H2 
luminosity as a function of the bulk kinetic luminosity of the Hi gas entrained 
in the flow. 

which is 10 3 — 10 4 times higher than the standard Galactic 
rate (Ch = 2.5 x 10" 17 s" 1 , Williams et al. 1998). This would 
require a cosmic ray pressure 1-10 times higher than the ther- 
mal pressure of the warm H2 gas (« 10 7 K cm" 3 ). Although 
we cannot exclude it for some of our sources, in the most ex- 
treme cases (3C 305, 3C 459, PKS 1549-79) that require the 
highest ionization rates, it is unlikely that such a departure 
from the equipartition of energy can be supported on galactic 
scales. 

Ogleetal. (2010) and Nesvadba et al. (2010) argued that 
the dissipation of the kinetic energy of the radio jet is the 
most probable source to power the H2 emission in radio MO- 
HEGs. The fact that all the radio sources studied here have 
bright H2 emission reinforces the idea that the radio jet, which 
drives the Hi outflow, is also responsible for the shock exci- 
tation of the warm H2 gas. As in 3C 326 (Ogle et al. 2007; 
Nesvadba et al. 2010), all Hi-outflow radio MOHEGs in our 
sample have jet kinetic powers that exceed the observed H2 
luminosity (see Table 2) and the Hi turbulent and bulk kinetic 
luminosities (see Table 14). The jet provides a huge reservoir 
of mechanical energy that can be dissipated over timescales 
of the order of 10 8 yrs, perhaps exceeding the jet lifetime 
(Nesvadba et al. 2010). Since the cooling time of the warm H2 



gas (« 10 yr) is much shorter than the dynamical timescales 
of the injection of mechanical energy (w 10 7 yr), the large 
masses of warm H2 gas requires that the gas is repeatedly 
heated. Spitzer H2 observations thus imply an energy cascade 
converting bulk kinetic energy of the multi-phase ISM to tur- 
bulent motions of much smaller amplitude within molecular 
gas (Guillard et al. 2009). This is what motivated us to model 
the observed H2 line fluxes with shock models (sect. 4.1). 
The turbulent heating of the molecular gas in powerful AGN 
is also suggested by observations of highly excited CO lines 
(e.g. Papadopoulos et al. 2010). 

6.2. Efficiency of the transfer from the jet kinetic energy to 
the Hi and H2 gas 

Two main mechanisms may be responsible for driving the 
observed outflow of Hi and ionized gas: the radiation pres- 
sure from the AGN (or the starburst if present), or the me- 
chanical impact of the radio source. For the sources studied 
here, based on the spatial correspondence of the Hi broad ab- 
sorption, optical outflowing component, and an off-nucleus 
bright spot in the radio jet for some of the sources, we favour 
the second mechanism. Outflows are most likely driven 
by the interaction between the expanding radio jets and the 
ionized/atomic gaseous medium enshrouding the central re- 
gions (e.g. Emonts et al. 2005). This may also be the case 
in even more powerful, high-redshift quasars (Nesvadba et al. 
2011b; Polletta et al. 2011) and high-redshift radio galaxies 
(Humphrey et al. 2006). 

The details of the physics that govern the energy transfer 
from the radio jet to the molecular gas are not yet understood. 
Being collimated, the jet itself cannot affect a significant frac- 
tion of the volume of the ISM of the host galaxy. Instead, 
the expanding cocoon of hot and tenuous gas inflated by the 
jet (e.g. Begelman & Cioffi 1989) could transfer part of its ki- 
netic energy to the molecular gas by driving shocks into the 
dense clouds or/and by turbulent mixing between hot and cold 
gas, similarly to the process proposed by Guillard et al. (2009) 
to explain the powerful H2 emission from the Stephan's Quin- 
tet galaxy collision. It also possible that part of the bulk ki- 
netic energy of the outflowing Hi gas is transferred to the 
molecular phase through shocks. The top (respectively bot- 
tom) panel of Figure 8 compares the H2 luminosity to the tur- 
bulent (resp. bulk) Hi kinetic luminosity. In all sources, the 
bulk Hi kinetic luminosity exceeds by more than one order of 
magnitude the observed H2 luminosities. The dissipation of a 
small fraction (< 10%) of this kinetic energy in the molecular 
gas could power the H2 emission. 

Numerical simulations of the impact of jets on galaxy evo- 
lution stress the importance of considering the multi-phase 
nature of the host galaxy ISM (e.g. Sutherland & Bicknell 
2007). Wagner & Bicknell (201 1) showed that the efficiency 
of transfer of kinetic energy and momentum from the jet to 
the dense («h > 10 2 cm" 3 ) gas can be high (10-70%), with the 
momentum flux of the clouds exceeding that of the jet. How- 
ever, the mid-IR IRS spectroscopy of the H2 rotational lines 
presented in this paper indicates that most of the molecular 
gas does not share the kinematics of the outflowing atomic 
and ionized gas. If most of the H2 gas was in the outflow, 
we would have observed blue-shifted H2 lines, similar to the 
ionized gas lines. This suggests that the dynamical coupling 
between the molecular gas and the more tenuous outflow- 
ing gas is weak, probably because of the high density con- 
trast between these phases. For density contrast of 10 4 , it 
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takes «4x 10 yrs to accelerate a 10 pc cloud up to a ve- 
locity of 500 km s" 1 (e.g. Klein et al. 1994), which is likely 
to be longer than the jet lifetime. If some H2 gas is effec- 
tively outflowing, the wind would remove this gas gradually 
over an AGN lifetime. It is very unlikely that the entire mass 
of the molecular disk would be entrained in the wind. Thus 
the fraction of the molecular gas expelled from these galaxies 
could be small, which does not fit with the current assump- 
tions made in the cosmological galaxy evolution models that 
assumes that AGN feedback sweep up most of the gas during 
the epochs of strong star-formation (e.g. Hopkins et al. 2006; 
Narayanan et al. 2008). 

The current numerical simulations show a very efficient 
transfer of momentum essentially because the dense («h > 
lcm -3 , T = 10 4 K) gas cools very fast. These calculations do 
not include the turbulent cascade and dissipation of the kinetic 
energy due to the supersonic turbulence within the dense gas 
phase, which would make the gas cooling longer, and there- 
fore the dynamical coupling less efficient. The dynamical 
coupling would be much more efficient if a significant frac- 
tion of the warm H2 gas is formed in situ, i.e. in the out- 
flow, from shocked Hi gas. On the contrary, our data suggest 
that the bulk of the warm H2 mass is difficult to entrain in 
the outflow, and is rather perhaps settled in the galactic disk. 
A similar conclusion is reached based on near-IR IFU obser- 
vations of nearby Seyfert galaxies (Storchi-Bergmann et al. 
2009; Riffel & Storchi-Bergmann 201 1), where the H2 gas is 
concentrated in the galactic plane (with a velocity structure 
consistent with a rotating disk), whereas the ionized gas is 
mostly distributed along the outflowing cone. 

On the other hand, in some of the sources studied here, the 
H2 lines are broad, showing that the H2 gas is very turbulent. 
This suggests that the molecular gas is heated by the dissi- 
pation of supersonic turbulence, which can prevent it from 
being gravitationally bound (and thus from forming stars), 
rather than expelled from the galaxy. In some cases, these 
turbulent motions are unlikely to be rotation-supported. As- 
suming the galaxies are pressure-supported with an isother- 
mal mass profile, we derived the escape velocities, v esc , from 
the stellar masses: v esc = \J2G M+/(5 r e ). We assume an ef- 
fective radius of r e = 2kpc. We find escape velocities that 
range from 70 km s" 1 (OQ 208) to 330 km s" 1 (4C 12.50). In 
all of the sources, the observed H2 velocity dispersions (up to 
w 300 km s" 1 , broader than that observed in nearby Seyferts) 
are comparable to their escape velocities, except for IC 5063, 
where er(H2) «3x v esc . Thus, a fraction of the warm H2 gas 
could become unbound and escape from these galaxies be- 
cause of its high velocity dispersion. However, the fraction of 
the H2 S(l) line emission having velocities larger than the es- 
cape velocity is small (< 10%), so it is unlikely that this pro- 
cess would strongly affect and empty the reservoir of molec- 
ular gas in these galaxies. 

We are complementing the observations presented here 
with higher resolution spectroscopy of near-IR H2 lines, inter- 
ferometric CO line observations, and Herschel spectroscopy 
and photometry, to complete our census of the physical and 
dynamical state of the molecular gas in these sources. For 
extended objects, the use of near-IR IFUs is necessary. Hope- 
fully, JWST/MIRI will allow us to probe mid-IR H 2 lines at a 
much better sensitivity and 5 x better spectral resolution than 
the IRS. 

7. CONCLUSIONS AND FINAL REMARKS 



We present mid-IR Spitzer spectroscopy of 8 nearby radio- 
galaxies where fast outflows of ionized and atomic gas were 
detected by previous optical and radio observations. Our main 
results and conclusions are the following: 

• In all of the sources, we detect high equivalent width 
H2 line emission from warm (100-5000 K) molecular 
gas, implying warm H2 masses ranging from 10 8 to 
10 10 Mq. The observed H2-to-PAH luminosity ratios 
are all above the values predicted from photo-ionization 
models, and we suggest that the H2 emission is asso- 
ciated with the dissipation of a fraction of the kinetic 
energy provided by the radio jet. 

• In 5 sources (3C 236, 3C 293, 3C 459, 4C 12.50 and 
PKS 1549-79), we securely detect blue-shifted wings 
(up to 3000 km s" 1 ) on the [Nell] (and [Nelll] for 
the highest SNR spectra) that matches remarkably well 
with the blue-shifted, broad Hi absorption associated 
with the outflow. 

• All but one of these sources have resolved and 
very broad mid-IR rotational H2 lines with FWHM> 
500 km s -1 (and up to 900 km s" 1 ), compared to only 
2% overall amongst the 298 AGN with such spectra 
in the Spitzer archive (Dasyra & Combes 2011). This 
suggests that the efficiency of kinetic energy deposition 
into the molecular ISM is higher in radio jet sources 
than in other types of AGN. 

• The kinematics of the warm H2 gas do not follow 
that of the ionized or Hi gas. The rotational H2 lines 
are systemically narrower than the mid-IR ionized gas 
lines, and do not exhibit asymmetric profiles with blue- 
shifted wings (except perhaps tentative detections in 
three targets, 4C 12.50, 3C 459, and PKS 1549-79). We 
conclude that, although very turbulent, the bulk of the 
warm H2 mass is not entrained in the wind. 

• We show that UV, X-ray and cosmic ray heating are un- 
likely to be the dominant source of H2 excitation. We 
argue that the dissipation (via supersonic turbulence) of 
a small fraction (< 10%) of the mechanical energy pro- 
vided by the radio jet can power the observed emission. 



This work is based primarily on observations made with the 
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Figure 9. Spitzer IRS high-resolution spectra of the Hi-outflow radio galaxies, zoomed on the wavelength coverage of the SH (Short High) module. 
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Figure 10. Detected spectral emission lines. The grey area is the Spitzer IRS instrumental profile (a Gaussian with a FWHM corresponding to the resolution of 
the IRS), and the red dashed line is the result of the Gaussian fitting of the line. The [OlV, Fell] blend profile has been fitted with two Gaussian components when 
both lines are detected. For lines showing wings, the broad component of the Gaussian decomposition is overplotted (orange dashed-dotted line). 



22 



P. Guillard et al. 









■ [Ne II] 


A 1000 km/s 






[Ne ill] 


1000 km/s 




15.4 15.5 15.6 15.7 
k rest On] 

IC5063 


[Ne V] 






14.20 M. 25 14.30 14.35 14.40 14.45 
\ rest O") 

IC5063 


Ne V 


1000 km/s 









0.08 






— 0.04 


^Ne^^ 




c 0.02 
0.00 






-0.02 






\ rest [urn] 


OQ208 




OQ20S 


" H 2 s 0) A 1000 km/s 


0.06 






T 0.02 






-0.02 




16.9 17.0 17.1 17.2 




15.4 15.5 15.6 15.7 



PKS1549 
; S(1) » 10 00 km /s 



[0 IV, Fe llJ 10 00 km /s 





0.020 - H z s (- 3 ) 
0.015 ■ 




[Ne II] 


1000 km/s 




12.70 12.75 12,80 12.85 12.90 12.9 
* rest [fun] 

PKS1549 


[Ne III] 




15.4 15.5 15.6 15.7 
* rest [urn] 

PKS1549 


[Ne V] 


1000 km/s 


14.2014.25 14.3014.35 14.4014.45 

k rast O") 
PKS1549 


Ne V i 


1000 km/s 




[0 IV, Fe II] iooo_ 




9.60 9.65 9.70 9.75 



Figure 11. Detected spectral emission lines. The grey area is the Spitzer IRS instrumental profile (a Gaussian with a FWHM corresponding to the resolution of 
the IRS), and the red dashed line is the result of the Gaussian fitting of the line. The [OlV, Fell] blend profile has been fitted with two Gaussian components when 
both lines are detected. For lines showing wings, the broad component of the Gaussian decomposition is overplotted (orange dashed-dotted line). 
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Figure 12. Hi excitation diagrams with two or three temperature component fits overlaid. The H2 temperatures, column densities, and masses are listed in 
Table 5. 
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Figure 13. H2 

and H2 masses 



excitation diagrams fitted with a combination of two C-shock models, for a pre-shock density of «h = 10 3 cm 3 . The shock model parameters, 
are listed in Table 6. 
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Figure 14. H2 

and H2 masses 



excitation diagrams fitted with a combination of two C-shock models, for a pre-shock density of «h = 10 4 cm 3 - The shock model parameters, 
are listed in Table 7. 
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Figure 15. Gaussian decompositions of the Hi absorption profiles observed with the Westerbork telescope. The data is from Morganti et al. (2005b), except for 
3C 293, where we used more recent observations taken by B. Emonts with the Westerbork telescope. 



